POWER SERIES SOLUTIONS OF NON-LINEAR 
^-DIFFERENCE EQUATIONS AND THE 
NEWTON-PUISEUX POLYGON 



J. CANO AND P. FORTUNY AYUSO 

Abstract. Adapting the Newton-Puiseux Polygon process to nonlin- 
ear q— difference equations of any order and degree, we compute their 
power series solutions, study the properties of the set of exponents of 
the solutions and give a bound for their q— Gevrey order in terms of the 
order of the original equation. 



1. Introduction 

The Newton Polygon construction for solving equations in terms of power 
series and its generalization by Puiseux has been successfully used countless 
times both in the algebraic [23] , |24j |17j and in the differential contexts [15] , 
021 Ch. V], [IE], [H], [7], [8], [TO], [32] (this is just a biased and briefest of 
samples, see also [9] and [TTJ Sec. 29] for an interesting detailed historical 
narrative). We extend its use to q— difference equations. 

Although this construction is primarily intended to give a method for 
computing formal power series solutions, we will use it for proving the q- 
analog of some results concerning the nature of power series solutions of non 
linear differential equations. Namely, we show properties about the growth 
of the coefficients of a power series solution (Maillet's theorem) and about 
the set of exponents of a generalized power series solution. 

The method allows us, first of all, to show that the set of exponents of 
any generalized power series solution of a formal ^-difference equation is 
finitely generated as a semigroup (in particular, it has finite rational rank 
and if the exponents are all rational, then their denominators are bounded). 
This mirrors the results of D. Y. Grigoriev and M. Singer in pH] for differ- 
ential equations. When the ^-difference equation is of first order and first 
degree, we give a bound for this rational rank (see Theorem [3] for a precise 
statement). We also study properties related to what we call "finite deter- 
mination" (Definition [4]) of the coefficients of the solutions. This is one of 
the places in which the case 1^1 = 1 is essentially different from the general 
case. For 1^1 ^ 1, we prove the finite determination of the coefficients. 

Maillet's theorem [20] is a classical results about the growth of the co- 
efficients a>i of a formal power series solution of a (non-linear) differential 

Date: September 4, 2012. 

2010 Mathematics Subject Classification. 39A13. 

1 



NEWTON-PUISEUX POLYGON OF NON-LINEAR g— DIFFERENCE EQUATIONS 2 



equation: it states that \m\ < i\ s R l , for some constants R and s. Among 
the different proofs (for instance (2Ql EH IS] ) , B. Malgrange's [22] includes 
a precise bound for s. This bound is optimal except for one case: when 
the linearized operator along the solution has a regular singularity and the 
solution is a "non-regular solution" , for which any s > works (see the last 
remark in Malgrange's paper); we shall refer to it as the (RS-N) case. In [7J, 
the Newton Polygon method allows the author to prove Maillet's result and 
to show convergence (i.e. s = 0) in the (RS-N) case. 

The first studies on convergence of solutions of non-linear ^-difference 
equations are due to Bezivin [3], [3] and [5]. The q-analog of Maillet's theo- 
rem states that when \q\ > 1, a formal power series solution of a ^-difference 
equation with analytic coefficients is q-Gevrey of some order s (see Defi- 
nition [5]). Zhang [33j proves this adapting Malgrange's proof to the case 
of ^-difference-differential convergent equations. In this paper, the adapta- 
tion of the Newton Polygon to ^-difference equations allow us to give a new 
proof of the g-analogue of Maillet's theorem and to extend it to the g-Gevrey 
non-convergent case. The bounds obtained for convergent equations match 
Zhang's in general and are more accurate in the (RS-N) case. However, we 
cannot prove convergence in this case unlike for differential equations. 

Notice that the "Newton Polygon" construction used in the case of linear 
operators by Adams pQ, Ramis [26], Sauloy [30] and others is different from 
the one presented here. In the linear case, the Newton Polygon is used to 
find local invariants of the operator while our Newton Polygon is constructed 
with the aim of looking for formal power series solutions. In Section [4] we 
describe the relation between Adams' Newton Polygon and Zhang's bounds. 
Adams' construction is also used in |19] to give conditions for the conver- 
gence of the solution(s) of analytic nonlinear (/-difference equations. 

For the reader's convenience, we include a final section with a detailed 
working example describing most of the constructions and the evolution of 
the Newton Polygon as one computes the successive terms of a solution. 

2. The Newton-Puiseux Polygon process for ^-difference 

equations 

Let q be a nonzero complex number. For j £ Z, let us denote by a 1 the 
automorphism of the ring C[[x]] of formal power series in one variable given 
by a J (y{x)) = y(q J x), that is, 

00 00 
aC^aiX 1 ) = ^2q ij a i x i . 

i=0 i=0 

Let P(x, Yq, Yi, . . . , Y n ) G C[[x, Yq, . . . , Y n ]] be a formal power series. For 
y G C[[x]], with oid x (y) > 0, the expression P(x,y,a 1 (y), . . . ,a n (y)) is a 
well-defined element of C[[x}] that we will be denoted by P[y\. We associate 
to P(x, Yq, Yi, . . . , Y n ) the q— difference equation 

(1) P(x,y,a 1 (y),...,a n (y)) = 0. 
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We will look for solutions of equation Q as formal power series with real 
exponents. We restrict ourselves to the Hahn field C((x w )) of generalized 
power series, that is, formal power series of the form X^ 7 eR c i xl whose 
support {7 I c-y ^ 0} is a well-ordered subset of M and c 7 G C. Hahn fields 
were essentially introduced in [16]; see [29] for a detailed proof of the ring 
structure and [31J for a modern study in the context of functional equations. 
We fix a determination of the logarithm and extend the automorphism a to 
C((x R )) by setting 

For y G C((x K )), its order ord(y) is the minimum of its support if y 7^ 
and ord(0) = 00. In subsection |2.2[ we shall see that if ord(y) > then 
the expression P(x, y, . . . , a n (y)) is a well-defined element of C((x R )), 

hence equation makes sense in our setting. 

Although we look for solution in the Hahn field, their support has some 
finiteness properties, as in the case for differential equations. We say that 
y G C((x R )) is a grid-based series if there exists 70 G M and a finitely 
generated semigroup V C M> such that the support of y is contained in 
70 + r. Puiseux series are the particular case of grid-based series in which 
70 G Q and T C Q. Puiseux series and grid-based series form subfieds of the 
Hahn field denoted respectively by C((x^)) 9 and C((x K )) 9 . We have 

C[[x}} C C((x Q )) g C C((x R ))f C C((x M )). 

If equation ([l]) is algebraic, i.e. of the form P(x, y) = 0, then by Puiseux's 
Theorem all its formal power series solutions are of Puiseux type. This is no 
longer true if instead of C, the base field is of positive characteristic, as the 
following example (due essentially to Ostrowski) shows: the equation — y p + 
xy + x = over the field Z/pZ has as solution the generalized power series 
V = X/Si w it n Mi = (p l ~ ~P % )- Notice that the exponents 

are rational but they do not have a common denominator and moreover 
fix < fi2 < ■ ■ ■ < — 1) so that they do not even go to infinity. Hence y 
is neither a Puiseux series nor a grid-based series. 

As in the case of differential equations, the number of generalized power 
series solutions of a given equation ([T]) is not necessary finite, neither all 
of its solutios are of Puiseux type. For instance, the g-difference equation 
Yq Y2 — Y"] 2 = has c x 11 as solutions for any c G C and 

2.1. The Newton Polygon. Let 1Z = C[[x R ]] be the subring of general- 
ized power series with non-negative order. For a finitely generated semi- 
group of T C M>o, the ring C[[x r ]] formed by those generalized power series 
with support contained in T is denoted by 1Zr- Let P G TZ[[Yq,Yi, . . . , Y^]] 
be a nonzero formal power series in n + 1 variables over 1Z. For p = 
(fto, Pi, ■ ■ ■ , Pn) e N n+1 , we shall write Y? = Y Q p0 ■ Yf 1 ■ ■ ■ Y? n ; we shall 
also write K[[Y}} instead of H[[Y , Y 1} . . . , Y n }]. The coefficient of Y p in P 
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will be denoted P p (x) G TZ and, for a£R, the coefficient of x a in P p (x) will 
be denoted P a , p G C, so that 

p= E and w= E p ^ x °' 

where for each p, T p is a well-ordered subset of M>o- We associate to P its 
c/oud of points C(P): the set of points (a, j) G M 2 for which there exists 

p = (p , ■ ■ ■ , p n ) with j = \p\ = po + pi H h p„ and P Q;P / 0. 

The Newton Polygon M(P) of P is the convex hull of 

C(P) = {(a + r,j) | (a,j) G C(P), r G M> }. 

A supporting line L of C(P) is a line such that C(P) is contained in the 
closed right half-plane defined by L and L n C(P) is not empty, that is a line 
meeting M(P) on its border. 

It will be convenient to speak about the co-slope of a line as the opposite 
of the inverse of its slope, the co-slope of a vertical line being 0. In order 
to deal with the particular case in which P is a polynomial in the variables 
Yq, Yi, . . . , Y n we shall make use of the symbol p_i(P), denoting — oo if P 
is a polynomial and otherwise. Hence the statement "p > p_i(P)" means 
"either p > or P is a polynomial" . 

Lemma 1. Let P G 7£[[Y]]. For any p > p_i(P) there exists a unique 
supporting line ofC(P) with co-slope p and the Newton polygon N{P) has 
a finite number of sides with co-slope greater or equal than p. If P is a 
polynomial then M(P) has a finite number of sides and vertices. If P G 
7^r[[Y]] f or some finitely generated semigroup T C M>o, then the Newton 
Polygon M(P) has a finite number of sides with positive co-slope. 

The unique supporting line with co-slope p will be denoted henceforward 
L(P;p). 

Proof. If P is a polynomial, let h be its total degree the variables Yq, . . . , Y n . 
Otherwise we define h as follows: since P / the set C(P) is nonempty; take 
a point q G C(P) and let L be the line passing through q with co-slope p. 
Let (0, h) be the intersection of L with the OY-axis. For each p G N n+1 , 
write a p = ord P p (x). Only the finite number of points (a p , |p|) with |p| < h 
are relevant for the definition of the line L(P; p) and for the construction of 
sides with co-slope greater or equal than p of Af{P). This proves the two 
first statements, the last one is a consequence of the fact that for a given 
a > 0, the set T n {r < a} is finite. □ 

For p > p_i(P), define the following polynomial in the variable C: 

<WC) = E Pa^q^C^, 
{a,\ P \)eL{P;ri 
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where w(p) = pi + 2p2 + • • • + np n . For a vertex v of M(P), the indicial 
polinomial is 

(a,|p|)=D 

2.2. Composition. For so,...,s n G C((x R )), the expression P(so, • • • , s n ) 
can be given a precise meaning under certain conditions. We consider on 
C((x R )) the topology induced by the distance d(f,g) = exp(— ord(/ — g)) 
which is a complete topology. 

If P is a polynomial, P(so, . . . ,s n ) is well-defined because C((x R )) is a 
ring. Otherwise, we impose ord(sj) > 0, for all i. Let p = mino<j<„{ord(si)}. 
For M G N, consider the polynomial P<m = Y1\ p \<m Pp( x ) Y p - The se- 
quence P<m(sq, ■ ■ ■ j s n ), M £ N, is a Cauchy sequence because the or- 
der of P p (x)sq° • • • Sn" is greater than or equal to p\p\. Its limit is pre- 
cisely P(so, . . . , s n ). Notice that if P £ lZr{[Y]] and all Sj G TZr, then 
P(s , ■ ■ ■ ,s n ) £ TZ r - 

Given so, . . . , s n as above, we define the series 

(2) P(s +Y ,...,s n + Y n ):= - ] ^^---^n)y P , 

where p\ = p \---p n l and = gyP0 fn P aY Pn ■ For generalized power 

series sq, . . . , s n with positive order it is straightforward to prove that the 
evaluation of the right hand side of ^ at sq, . . . , s~ n is P(so + sq, . . . , s n + s n ). 

If y e C((x R )) has ord(y) > P-i(P), then P(y,a(y), . . . ,a n (y)) is well 
defined because ord(cr fc (y)) = ord(y). We also remark that if y G 7^r) then 
(7 k (y) G TZr- The following notations will be used in the rest of the paper: 

P[y]=P(y,a(y),...,a n (y)), 

(3) P[y + Y] =P(y + Y ,a(y)+Y 1 ,...,a n (y)+Y n ). 

We are also going to make use of the little-o notation: o{x p ) will mean a 
generalized formal power series with order greater than p or the zero series 
if p = oo. The following is essentially what motivates the Newton polygon 
construction: 

Lemma 2. Let y = ci* 1 + o(x At ) G C((x R )), and p > p-\(P). Let (V, 0) be 
the intersection point of L(P; p) with the OX -axis. Then 

P[y} = $ (P ., tM) (c)x» + o(x»), 

In particular, if y is a solution of the q-difference equation |Ip then 

$(P; M )( C ) =°- 

Proof. We may assume that P = P<m for some M, otherwise let M G N 
such that Mp > v. The truncation of P[y] up to order v is equal to that of 
P< M [y] and also $(p ;m )(C) = $( P < m;m) (C). 
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Write a p = ord P p , and notice that v is the minimum of a p + p \p\, for 
p 6 N n . The following chain of equalities proves the result 

P<m[cx" + o(/)] = 

E P P {x) (c/ + o(/)r(g fl c/ + o(x' 1 )f ■•■(^cZ + ot/))''" = 

|p|<M 

2 { p w x " P + °(^ p )} {c^g^^x^ 1 +o(x^)} = 

|p|<M 

E {Pa p ,pcW tf"^) X <*P+P\P\ +o(x a P + ^)\ = 
\p\<M 

E p «-p C ' P| | xU + °0O = *(P;ri( c ) + °(^)- 

The last equality holds because L(P; p) = {(a, b) \ a + fib = u}. □ 

Let y S C((x R )) be a generalized power series and S be its support. If S 
is finite, denote by io{y) the cardinal of S, otherwise oj(y) = oo. Consider 
the sequence fii € S defined inductively as follows: po is the minimum of 
S and for < i < co(y), pi + \ is the minimum of S \ {po,pi, . . . ,Pi}. Let 
Ci G C be the coefficient of in y. 

Definition 1. We shall call the first oj terms of y to the generalized power 
series £o<i< u (v) d x w . 

Notice that if the support of y is finite or has no accumulation points then 
y coincides with its first oj terms. 

Corollary 1. Let y be a solution of the q-difference equation |Ip and let 
'YltiCiX^ be the first oj terms of y. Let Pi be the series defined as: 

P := P, and P t+1 := P[c z x^ + Y], 0<i< u{y). 

Then, for all < i < co(y) , one has 

®{P^){ c i) = 0, and Mi-l < Mi) 
where we denote [i-i = p-\(P). 

Proof. Let y~k = U — Y2i=o c * xA ^> then Pk[y~k] = and the first term of y^ 
is Cfc x^ k . □ 

Let P G 7£r[[y]] an d let £°^ c$ x^ be a series with /i_i(P) < m < fii+i, 
for all < i < oo (We do not impose that Cj 7^ 0, but the sequence is 
strictly increasing). Consider the series Pq := P and -Pj+i := Pi[ciX^ + Y]. 

Definition 2. We say that £°^ c * x ^ sa tisfies the necessary initial condi- 
tions for P , in short NIC(P), if $(p j;w ) (c*) = 0, for all i > 0. 
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The above Corollary states that the first u terms of a solution of P[y) = 
satisfy NIC(P). In this section and the next one we shall prove in Propo- 
sition [2] the reciprocal statement for P G T^rll^]]: if ]C£o c i x ^ satisfies 
NIC(P), then lim^oo /ij = 00 and c i x ^ ^ s an actual solution of the 

(/-difference equation P[y] = 0. This implies in particular that solutions of 
P[y] = coincide with their first oj terms. 

A method for computing all the series satisfying NIC(P) with a ^ 0, for 
all i, is the following one: 

Procedure 1 (Computation of a power series satisfying NIC(P)). 
Set Pq := P and pL-\ := /i_i(P). 

For i = 0, 1, 2, . . . do either (a.l) or (a. 2) and (b), where: 



(a. I). If y = is a solution of Pi[y] = 0, then return X)fc=c 



k=0 C k x 

(a. 2). Choose /i, > [Ai-i, and 7^ Cj G C satisfying $(p. ^ (cj) = 0. 
If neither (a.l) nor (a. 2) can be performed then return fail, 
(b). SetP l+1 (Y) := P[ Ci x»> +Y]. ' 

If fail is returned at step k of the above Procedure, this means that 
there are no a solutions of P[y] = having X^fco c i xfli as its ^ rs t ^ terms. 
To prove this, assume that z is a solution having X^fco c i x>li as its first fc 
terms. Either z = Y2i=o c i x>ii 1 i n which case y = would be a solution of 
Pk[y] = and (a.l) would have been performed, or z — Yli=o c i x>li would 
have a first term of the form CkX^ k so that (a. 2) could have been performed. 

In order to carry out (a. 2) in the above Procedure, one has to deal with 
the following formula with quantifiers 

(4) 3/i>//, 3ceC,c/0, $ (P;a0 (c) = 0. 

The Newton Polygon provides a way to eliminate the quantifiers. Fix // > 
//_l(P); by Lemma [TJ N(P) has only a finite number of sides L±, L2, ■ ■ ■ , Lt 
with co-slopes greater than fi' . Let 71 < 72 < . . . < "ft be their respective 
co-slopes and denote by Vi-i and t>j the endpoints of Li. Take fi > fi' . 
Either fj, = jj for some 1 < j < t, or jj < < 7j+i for some < j < i 
(writing 70 = // and jt+i = 00). If /i = jj, then L(P;yu) n M(P) = Lj 
and $(p ;/Lt )(C) depends only on the coefficients P ajP of P with (a, \p\) G Lj. 
Otherwise, jj < fi < Jj+i for some j and L(P; /i) n AT(P) is just the vertex 
■Uj = (a, 6), so that <3?(p ;At )(C) = C b ■ ^ (p.„ . Thus, there exists c / 
and with 7, < /x < 7^+1 such that <3?(p ;At )(c) = if and only if there exists 
p, with jj < fx < 7j + i and ty(p. v ^ {q^) = 0. This proves that equation (Ij) is 
equivalent to the quantifier- free formula obtained by the disjunction of the 
following formulae: 

(5) * (P . 7 . ) (c) = 0, l<j<t, 

(6) *(p ;% )( T ) = °> M = log T/ log q, 7i < M < 7j+i, < j < t. 

2.3. The pivot point. For P G T^r-flX]] an d > /z_i(P), we shall denote 
by Q{P;n) the point with highest ordinate in L(P; /i) fl N{P). For P = 
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P[cx M + Y] (as in equation (|3])), the following Lemma describes the Newton 
Polygon of P: 

Lemma 3. Let h be the ordinate of Q(P; p) and consider the half-planes 
h+ = {( a ,b) G M 2 | b > h}, hr = {(a, b) G R 2 \ b < h}. If L(P;p)+ is 
the closed right half plane defined by L(P; p) and (y, 0) is the intersection 
of L(P; p) with the OX -axis, then 

(1) Af(P) n h + = Af(P) n h + , in particular Q(P; p) G Af(P). Moreover, 
for any a and p with (a, \p\) = Q(P;p), the coefficients P a ,p o,nd 
Pa,p are equal. 

(2) M(P) n hr c L(P ; p)+ n hr, 

(3) The point (v, 0) G N(P) if and only if $(p ;(U )(c) / 0. 

Proof. Write M p (Y) = P p (x)Y p and a p = oidP p (x). It is straightforward 
to show that M p [cx^ + Y]= M p (Y) + V(Y) for some V(Y), whose cloud 
of points is contained in the set A p = {(a, b) \ b < \p\} Pi L{M p \ p) + . This 
proves part (2). If Q = (a, p) belongs to M{P) n h + , then there are no 
points Q' = (a' , p') G M(P), except Q itself, such that Q G >Ly. This proves 
part (1). Part (3) is a consequence of Lemma [2j □ 

Corollary 2. Let p, > p. Then either Q{P;p) = Q(P,p) or the ordinate 
ofQ(P,p) is less than the ordinate of Q(P; p). If <J?(p ;M )(c) ^ 0, then the 
ordinate of Q(P;p) is zero. 

Proof. The previous Lemma implies that Q(P; p) is a vertex of M{P) and 
L(P;p) = L(P;p). Hence Q(P;p) = Q(P;p). Since p > p, Q(P;p) is a 
vertex with ordinate less than or equal to the ordinate of Q(P; p) = Q(P\ p)- 
For the second part, assume that <I?(p ;A1 )(c) 7^ 0. By the same Lemma, the 
point (^,0) G M{P), so that the segment whose endpoints are (^,0) and 
Q(P; p) is the only side of M{P) with co-slope greater than or equal to p, 
from which follows that Q(P; p) = (^,0). □ 

Let P G Kr[[Y}] and take a series tp(x) = Ei^o c * x ^ with M-i(-P) < 
Pi < pi + \ for all < i < 00. (Notice that we do not impose that q 7^ 0, 
but the sequence (pi)ieN must be strictly increasing). Writing Pq := P and 
Pi + i := Pi[ciX Pi + y], let = Q(P{;pi). By the previous Corollary, the 
ordinate of Qi is less than or equal to the ordinate of Qi+i- Since these are 
natural numbers, there exists N such that for % > N, the ordinate of Qi is 
equal to the ordinate of Qn (it stabilizes). By the same Corollary, we know 
that actually Qn = Qi, for all i > N. This leads to the following 

Definition 3. The pivot point of P with respect to ip(x) is the point Q at 
which the sequence Qi stabilizes and is denoted by Q(P;ip(x)). We say that 
it is reached at step N if Qn = Q(P', 4>{x)). 

Let Qn = {ot, h) be the pivot point just defined. From part (1) of Lemma|| 
follows that {Pn)o,p = (Pi)a,p for all i > N, and for all p with |p| = h. In 
particular, the indicial polynomials ^(p-.q^^T) are the same for all i > N. 
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We shall say that the monomial Y p (resp. the variable Yj) appears effectively 
in the pivot point if (Pj\r)a,p 7^ (resp. for some p with pj > 0). 

Lemma 4. Let P and if)(x) = Y^=o Ci xMi ^ e as a ^ ove - The following state- 
ments are equivalent: 

(1) The ordinate of the pivot point of P with respect to X^So Ci xfli * s 
greater than or equal to one. 

(2) The series YT=o c i x ^ satisfies NIC(P). 

In case \i~mpi = oo, these statements are equivalent to 

(3) The series ip(x) is a solution of P[y] = 0. 

Proof. Assume statement (1). The ordinate of Qi+i is non-zero and by the 
above Corollary, $fp.. jU .)(cj) = 0, which proves (2). Assume now that state- 
ment (1) is false, so that the ordinate of the pivot point is zero. This means 
that there exists some N such that Qn has ordinate zero. By definition of 
Qn we have that L(P/v; p^) nJV(Pjv) is just the point Qn = (a,0). Then 
^(P N -,fi N )(C) 1S a non-zero constant (namely the coefficient of x a in P/v), 
therefore it has no roots, in contradiction with $(p n . /J/N \(cn) = 0. This 
proves the equivalence between (1) and (2). By CorollaryfTJ (3) implies (2). 

Assume (1) holds and that lim//, = oo. Write tpk(x) = Yli=o c i xfli an d 
notice that Pj = P[ipi(x) + Y], in particular, P[ipi(x)] = Pj[0] = (Pj)g- Let 
Q = (a,h) be the pivot point of P with respect to tp(x). Since L(Pf, pn) 
contains the point Q, ord(Pj)g > a + h/M and since h > 1, the sequence 
ord P[4>i(x)] tends to infinity and we are done. □ 

Corollary 3. Let Yli^o °i x>li be the first oj-terms of a solution of P[y] = 0. 
Then the pivot point of P with respect to J2i^o c * x ^ has ordinate greater 
than or equal to one. 

2.4. Relative pivot points. The above construction of the pivot point can 
be made relative to any of the variables Yj, < j < n, and more generaly, 
relative to any monomial Y r ', with r = (ro, n, • • • , r n ) G N™ +1 , as follows: 

Fix r G N™ +1 . The cloud of points of P relative to Y r is defined as the 
set C r (P) = {(a,\p\) | 3p, with P a p ^ 0, and r ■< p}, where r < p means 
that n < pi, for all < i < n. It is obvious that C r (P) C C(P). 

Assume that C r (P) is not the empty set, then we may define the line 
L r (P; p) as the leftmost line with co-slope [i having nonempty intersection 
with C r (P). The point Q r (P;p) will be the one with greatest ordinate in 
L r (P;p)nC r (P). 

If H denotes H = d Qyr , the cloud C r (P) is not the empty set if and 
only if H is not the zero series. In this case, consider the translation map 
r(a, b) = (a, b — |r|). It is straightforward to prove that C r (P) = r _1 (C(P)). 
Hence, L r (P; p) = r- 1 (L(P; p)), and Q r (P; p) = t- 1 (Q{H; p)). 

Let tp(x) = X^So c * x ^' w ith po > p-i(P), and assume the notation of 
Definition |3j Denote also Hq = H and Pj+i = H^CiX^ + Y\. By the chain 



NEWTON-PUISEUX POLYGON OF NON-LINEAR g— DIFFERENCE EQUATIONS 10 



rule, 

r) h P 

(7) W^ = Hl > 

The sequence of points QriPi',^) = T~ 1 (Q(Hf, /ij)) for i > stabilizes at 
some point denoted Q r (P; i^(x)) and which we call the pivot point of P with 
respect to ip(x) relative to Y r . Therefore 

(8) Q(H;i;(x)) = T(Q r (P;i;(x))) 

Remark 1. Since H / 0, then Hi ^ 0, for i > so that C r (Pi) is not empty, 
for i > 0. This proves that Q r (P; m) and Q r {P;i^{x)) are well-defined 
provided the monomial Y r appears effectively in P. 

From now on, we shall denote ej the vector (0, . . . , 0, 1, 0, . . . , 0) where the 
1 appears at position j + 1, for j = 0, . . . , n. Thus, ej = (<%)o<i<n £ N n+1 
where 5{j is the Kronecker symbol. 

Proposition 1. Let Q = (a, h) be the pivot point of P with respect to ip(x). 
Assume that the monomial Y r> appears effectively in Q. Let r € N n+1 ; with 
r ■< r' , and H = • Then the pivot point of H with respect to if)(x) is 
(a,h — \r\). In particular, if r = r' — ej, for some i such that r\ > 1, then 
the ordinate of the pivot point Q(H; ip(x)) is one. However, for r = r' , one 
has Q(H;ip(x)) = (a,0) and therefore tp(x) is not a solution of H[y] = 0. 

Proof. Assume the pivot point Q is reached at step N, thus Q G C r >{Pi) C 
C r (Pi) for all % > N. From C r (Pi) C C(Pi) and the fact that Q = Q(Pj; 
for all i > N, one infers Q = Q r (Pi',fJ.i) = Q r '(Pf, Hi) for all i > N. This 
means that Q is the pivot point of P with respect to ip{x) relative to Y r and 
also relative to Y r . As we have seen before, r r (Q) = (a,h — |r|) is the pivot 
point of H with respect to ip(x). The third statement is a consequence of 
Lemma 01 □ 

Corollary 4. Let ip(x) = Yl^o c * xfli ^ e a solution of P[y] = with lim /ij = 
oo. If the pivot point (P;ip(x)) has ordinate greater than one, then there 

exists a non trivial derivative H = ^yr- of P, such that ip( x ) ^ s a solution 
ofH[y}=0. 

Proof. Let Y r ' be a monomial that appears effectively in the pivot point 
Q = Q(P;tp(x)). Since Q has ordinate greater that one, r' can be chosen 
with \r'\ > 2. Let r be such that r < r' and 1 < \r\ < \r'\. By the 
Proposition, the pivot point of H with respect to ip(x) has ordinate greater 
than or equal to one. By Lemma |4j tfi(x) is a solution of H[y] =0. □ 

Lemma 5. Let Q{P;ip{x)) = (a,b) and Q r {P;^{x)) = (a',b') be respectively 
the general pivot point of ip{x) and the pivot point of ip{x) relative to Y T . 
If the sequence ^ of exponents of ip(x) tends to infinity, then b' > b and if 
b' = b, then a' > a. 
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Proof. Assume that both pivot points have been reached at step N. For 
any i> N, the point (a', b') belongs to the closed right half plane L(P; fii) + 
because C r (Pj) C C(Pj). Since (a, b) G L(P;/ij) for all i> N, and lim/ij = 
oo, the intersection of all the half planes L(P; /Uj) + for i > N, is the region 
formed by the points in L(P; hn) + with ordinate greater than or equal to b. 
This proves the Lemma. □ 



Throughout this section, we assume that T is a finitely generated semi- 
group of R>o and that P is a nonzero element of 72-r[[^]]- We also as- 
sume that q ^ 1: the case q = 1 is reduced to the case n = considering 
P(Yq, Yq, . . . ,Yq). This section is devoted to proving the following results: 

Theorem 1. If y G C((x K )) is a solution of equation then it is a grid- 
based formal power series. 

Proposition 2. If ip(x) = Y^o c i x ^ satisfies NIC(P), then tfj(x) is a 
solution of P[y] = 0. 

Definition 4. Let y £ C((x R )) and P £ ft r [[Y]]. We say that y is finitely 
determined by P if there exist positive integers k and h, such that if y^ de- 
notes the first k terms of y then y is the only element z 6 C((x K )) satisfying 
the following property: "z^ = y^ and Q[y] = if and only if Q[z] = 0, for 
any Q = 8 d yf i with \r\ < h." 

Theorem 2. If \q\ ^ I, then any solution y of equation |7]j is finitely 
determined by P. 



The hypothesis \q\ ^ 1 is necessary: let P = Yq — Y% and q = \J — 1. Any 

series Yli^o c 4* (^ or arbitrary constants c^) is a solution of P[y] = 0. 
Since dP/dYo[y] = and dP/dY\[y\ = have no solutions, and higher order 
derivatives of P are zero, none of these solutions is finitely determined by P. 

To prove Theorem [T] we shall proceed as follows: let tp(x) be the first u 
terms of y. By Corollary [T] above, ip(x) satisfies NIC(P). We reduce the 
equation to what we call quasi solved form using Lemma [7] below and in this 
case, by Lemma[8j we infer that supp ip(x) is contained in a finitely generated 
semigroup of M>o- By virtue of the following Remark [2j y coincides with 
ip(x). As a byproduct, we shall obtain a recursive formula for the coefficients 
of a solution which will be useful in the subsequent section. 

Remark 2. Let T be a finitely generated semigroup of M>o- For any real 
number k, the set m{r | r < k} is finite. Hence T is a well-ordered set with 
no accumulation points and its elements can be enumerated in increasing 
order: V = {7i}i>o, with 7; < 7*4.1 and liniTi = 00. Let ip(x) = J^iZo ^ 
be the first oj terms of an element y G C[[x R ]]. If supp ip(x) is contained in V 
then either it is finite or lim//j = 00. In both cases, y = ip(x). In particular, 
any element of C[[x R ]] whose support is contained in T coincides with its 
first co terms. 



3. FlNITENESS PROPERTIES 
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3.1. Quasi solved form. We say that the equation 
(9) P[y] = 0, ord(y) > 0, 

is in quasi-solved form if the point (0, 1) is a vertex of A/"(P) and (0, 0) G" 
C(P). If this is the case, let \&(T) be the indicial polynomial of P at (0, 1), 
S = {/i £ 1 *(^) = 0} and FJ+ = £ n M> . We say that equation @ is 
in solved form if S + is the empty set. 

Remark 3. The polynomial ^(T) can be written \I>(T) = Po,e + ^b,e a 7" + 
• • -+Po,e n P n G C[T]. Its degree m is the largest index such that the variable 
Y m appears effectively in the point (0, 1). If the equation is in quasi-solved 
form, vl/(T) is a nonzero polynomial because (0, 1) G C(P). If \q\ ^ 1, then 
S is finite. If case |(;| = 1 (and q ^ 1), then £ is the finite union of the sets 

^ r = arg(g) arg(g) ^' ^ or those complex roots r of ^>(T) with modulus one. 
Recall that we have fixed a determination of the logarithm to compute q^, 
hence arg(g) is also fixed. The following Lemma implies that S r n M>o is 
a finitely generated semigroup. Therefore S + generates a finitely generated 
semigroup of M>o- 

Lemma 6. Let 7 G M. and 71, 72, •• -,7s positive real numbers. Then the 
semigroup T o/M>o generated by the set A = (7 + 71 N + • • • + 7 S N) n M>o 
is finitely generated. 

Proof. Let A be the set of (m, . . . , n s ) G N s such that 7 + ]T) > 0. By 
Dickson's lemma, the number of minimal elements in A with respect the 
product order are finite. Hence T is generated by 71, . . . , 7 S and the family 
7 + n Hi > f° r au minimal element (m, . . . ,n s ) of A. □ 

3.1.1. Change of variable z = x^y. Let P G T^rfl^]] and 7 > //_i(P). Define 
P[x 7 y] as the series 

qi^p) X 1\P\ p p (x) Y^Yf 1 ■ ■ ■ Y? n G C((x R )Y [[Y]]. 

p 

If (u, 0) is the intersection point of L(P; 7) with the OA-axis, then all the 
coefficients of the series P[x 7 y] have order greater than or equal to v. De- 
fine 7 P = x~ l 'P[x J Y]. The coefficients of 7 P are in TZr*, where T* is the 
semigroup of M>o generated by (— v + V + 7N) n M>o- By Lemma [6j T* is a 
finitely generated semigroup of M>o- 

The transformation P 1— >• 7 P corresponds to the change of variable 2 = x^y 
in the following sense: for a series y, with ordy > 7 + fi-i(P), one has 
7 P[x~ 7 y] = x _1/ P[y], in particular, P[y] = if and only if 7 P[x~ 7 y] = 0. 

Let f (a, b) be the plane affine map f(a,b) = (a — v + 76, b) , which satisfies 
f(Cj(P)) = CjpP) for < j < n. In particular, f(N(P)) = N{~P), and 
f maps vertices to vertices and sides of co-slope [i > 7 to sides of co-slope 
fj, — 7. Moreover, f(L(P;^)) = L( 7 P; — 7), in particular f(L(P;j)) = 
L( 7 P; 0) is the vertical axis. Therefore, Q(P; /x) and Q( 7 P; /x — 7) have the 
same ordinate. Let X^So c « x ^ an d Pi be as in the definition of pivot point 
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(Definition [3]) . Assume 7 < /io and set H = 7 P, Hq = H and Hj+i = 
Hi[ciX^ i ~' 1 + Y]. It is straightforward to prove that 7 Pi = Hi, so that 
f(Q(Pj; /Xi)) = Q(Hi] fa — 7) and, in particular, they have the same ordinate. 
Then the image by f of the pivot point of P with respect to c i x/ii is 

the pivot point of 7 P with respect to Xl£o CixMi_7 anc ^ the same holds for 
relative pivot points. By Lemma [4j this implies that c i x ^ satisfies 

NIC(P) if and only if c^" 7 satisfies NIC( 7 P). 

Finally, if v G C(P) then f\v) G C( 7 P) and ^ {jp . f{v)) (T) = ty (P . v) (q"<T). 

Lemma 7. Assume that t/j(x) = c « satisfies NIC(P). Then there 

exist a finitely generated semigroup F* , a series P* G IZr* [[^o 5 Y\, ■ ■ ■ , Y n ]], 
an index N and a rational number 7 with fasr-i < 7 < /A/v> swc/i £/tai i/ie 

(10) P*[z] = 0, ordz>0 

is in quasi solved form and ip*{x) = Ya^Ln Ci x ^~^ satisfies NIC(P*). 

Proof. We may assume that the ordinate of the pivot point of ip{x) with 
respect to P is one. Otherwise, by Proposition [TJ we may replace P by any 
of its derivatives 7^7- , where the monomial Yj Y r appears effectively in the 
pivot point, for some j. We remark that the coefficients of any derivative of 
P also belong to IZr- Let Q = (a, 1) be the pivot point of P with respect to 
ip(x) and use the notation of Definition [3j Pq = P, P%+i = P^CiX^ + Y] and 
so on. In particular, let the pivot point be reached at step N' — 1 for some 
N', Consider any integer N > N' . Denote r = T and F i+1 = Fi + faN. 
Notice that the coefficients of P, belong to T^r, • 

Let 7 be a rational number such that fas-\ < 7 < fav and set P* = 
7 P/v G ^-r^[[^]]- Since the pivot point Q has been reached at step N — 1, 
Q G L(P N -r, pLN-i) nL(P N ;fi N ). By Proposition |3j Q G L(P N ; yUjv-i). 
Hence Q G £(P/\r; fas-i) H L(Pn; fay); since /xj\r-i < 7 < /Mr, we conclude 



that Q(P/v;7) = Q = (a, I). So, by the properties described in 3.1.1 the 
point (0, 1) is in C(P*), the equation P*[y] = is in quasi solved form and 
the pivot point of P* with respect tp*(x) is (0,1). By Lemma [4j ip*(x) 
satisfies NIC(P*). □ 



Lemma 8. Assume equation (10) is in quasi-solved form and let £(cc) = 
c i x>l1 > w ith Mo > 0, be a series satisfying NIC(P*). Then the support 
of £(x) is contained in the finitely generated semigroup F' = F* + S + N. In 
particular, either the support of £(x) is finite or lim fa = 00 and in both 
cases £(cc) is a solution of equation uT 



Proof. Let P = P* and P+i = Pfaxf* + Y] for i > 0. We first prove 
that Q{Pi\ fa) = (0, 1) for all i > 0. We do this showing, by induction on 
i, that jV(Pj) is contained into the first quadrant of the plane and that the 
point (0,1) G C(p). This holds for Pq because of the hypotheses on P*. 
Assume that the statement holds for Pj. Since fa > 0, the line L(Pi; fa) 
either contains the point (0,1), and then Q(Pf,fa) = (0,1), or L(Pf, fa) 
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meets Af(Pi) at a single point with zero ordinate which is Q(Pf, m). If the 
latter happens, from Corollary [2j we infer that the pivot point of P* with 
respect to £(x) has zero ordinate, in contradiction with the fact that 
satisfies NIC(P*). Hence Q(Pi] /Uj) = (0, 1). By Lemma[3j (0, 1) is a vertex of 
A/"(Pi+i) and since P{+i G 7£[[Y]], its Newton polygon is contained in the first 
quadrant. This proves the induction step and that Q(Pi] fa) = (0, 1), i > 0. 

The fact that Q{Pi] fa) = (0, 1) implies that the polynomial $(p. ;At .)(C) is 
equal to ^(q^^C + Coeff x Mi Y-), where ^(T) is the indicial polynomial 
of P at (0, 1) and Coeff^; Y^-) is the coefficient of x^ of Y§Y? ■■■Y® 
in Pi. Since <3>(p i;/J .)(cj) = because (,(x) satisfies NIC(P*), the following 
equations hold: 

(11) ^(q^)ci + Coeff(P i; x^y fi ) = 0, i > 0. 

Let us prove, by induction, that Pi G 7^r'[[^]]) for all i > 0, and that the 
support of is contained in T' . By hypothesis, Pq G TZr'[[Y]]. Assume 
that Pi G Ti TI [[Y}}. \ia = 0, then P i+1 = P { G K T > [[Y]] and fa supp(£(x)). 
If, on the contrary, Cj ^ 0, we can prove by contradiction that fa G V: 
assume that fa T', in particular /ij S + , hence ^(g^) 7^ 0. From 



equation (11), Coeft (Pi] x^ Y&) + 0, we deduce that fa G supp((Pj)o) Q V. 
So Pi + \ = PilciX^ + Y] belongs to T^r'IM] which proves the induction step. 

The set supp^(x) has no accumulation points in K because T' is a finitely 
generated semigroup of R>o and supp£(x) C T' and we are done. □ 



Corollary 5. Let y be a solution of equation (10) which is in quasi solved 

form. Let V = {7i}^ > 7* < 7*+i f or a ^ ^ en 2/ = ^ x7i w ^ 

satisfying the following recurrent formula: 

(12) f(g 7 'K = -Coeff(P*[(iia; 71 + ••• + d^x^}; x^) , i > 1. 



7/ S + is finite and z is another solution of equation (10) with ord(y — z) 
greater than any element of T, + , then y = z. 

Proof. Let £(x) be the first u terms of y. Then supp£(x) C r', and by 
Remark [2j y = £(x) G T^r'- Hence we may write y = diX^ because 

70 = and ordy > 0. Since £(x) satisfies NIC(P*), the same reasoning 
as in Lemma [8] up to equation (11) holds. The coefficient Coeff (Pf, x^'Y-) 



is equal to the coefficient of x" H in P*[d\x~ tl + • • • + di-ix^- 1 ], which gives 
equation ( p"2j ). To prove the last statement, write z = X^i^i x7i - ^ 7« 
is greater than any element of S + , then ^(q^) 7^ 0, and di is completely 
determined by d\, . . . , dj-i, so that y = z. □ 

Proof of Proposition \^ Applying Lemma[7]to ^(a?) we obtain equation (10), 
and applying Lemmaj8ko £(x) = YIULn Cj 3 ^" 7 we conclude that fa — 7 G V , 
for i > N. Since 7 > fa), the set (7 — fa) + L' is included in M>o- Let F" 
be the semigroup generated by (7 — fa) + V . By Lemma [6^ r" is finitely 
generated. Let V" be the finitely generated semigroup r"+^ i= Q (fa— fa)) N. 
The set suppV'(^) is contained in fa + T'", so that lim/ij = 00. By Lemma|4j 
V'(x) is a solution of P[y] = 0. □ 
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Proof of Theorem^ Let ip(x) = Y2i^o c i xfli be the first u) terms of y. By 
Corollary [TJ if)(x) satisfies NIC(P). As in the proof of Proposition [2j there 
exists a finitely generated semigroup V such that supp^(x) is contained in 
fio + T. By Remark [2j y = ip(x), so that y is grid-based. □ 

Proof of Theorem plj Let y be a solution of equation ([!]) . By Theorem [TJ y 
coincides with its first ui terms. Write y = YliLo CiX^ and let Q = (a, h) be 
the pivot point of P with respect to y. Apply Lemmas [7] and [8] to y: let N 
and 7 be as in Lemma [7j we may assume that the pivot point Q is reached 
at step N — 1. Since \q\ 7^ 1, S + is finite by Remark [3J Since lim/Uj = 00, 
there is k > N such that fx^ — 7 is greater than any element of E + . 

Consider z S C((x K )) with the same first k terms as y and satisfying that 
for any H = 7^77-, with \r\ < h, H[y] = if and only if H[z] = 0. We have 
to show that y = z. 

Since P[y] = 0, then P[z] = 0, and z coincides with its first uj terms. 
Write z = S^o^* x<5i - By hypothesis, = d{ and /Xj = 8{ for < i < k. 
Denote P^ = P, P' i+l = P^x 6 * + Y] and P = P and P i+1 = Pi[ciX^ + Y], 
for i > 0. Obviously, P{ = P!, for < i < k. In particular Q = Q{P^; = 

Q(P' N ;S N ). 

If the pivot point of P with respect to z is also Q, then apply Lemmas [7] 
and [8 to z in the same way as to y: choose the same derivative ^gy/" , 
the same iV and the same 7 to obtain the same P*. This can be done 
because Pj = P/, for < i < k. This implies that £(x) = X^^jv c i x>Mi ~' y an d 
C(^) = E t =7v diE 5 ' -7 both satisfy NIC(P*). By Corollary g £(x) = £(x), 
which implies y = z. 

Let us show by contradiction that the pivot point Q' of P with respect to 
z must be Q. Assume Q' / Q. Hence Q = Q{P' N ; 5m) is not the stabilization 
point of the sequence Q(P-;5i). This implies that Q' has ordinate h! < h. 
Let Y r be a monomial that appears effectively in the pivot point of P with 

respect to z, so that \r\ = b! . Let H = 9 dY r ' • By Proposition 1, 7^ 0; in 
particular H 0. We claim that JET[y] = 0. By Remark [TJ the pivot point 
Q r of P with respect to y relative to Y r is well-defined. Since lim /ij = 00, 
by Lemma [5J the ordinate of Q r is /i" > /i. The pivot point of H with 
respect to y has ordinate h" — h! > h — h! > 1 . By Lemma [TJ y satisfies 
NIC(P) and so H[y] = 0, which proves our claim and finishes the proof the 
the Theorem. □ 

3.2. Bounding the rational rank in the case of order and degree 
one. Recall that the rational rank of a semigroup S C K is the dimension of 
<S>, the Q- vectorial subspace of K generated by S. It is denoted rat. rk(5). 
In what follows T denotes a finitely generated semigroup of 7£>o, as above. 

Theorem 3. Assume \q\ 7^ 1. Let P = A(Yq)-\-B(Yq)Yi be a nonzero series, 
where A, B G 7^r[[^o]]- Let y be a solution of P[y] = 0, with ordy > /i_i(P). 
Then rat. rk(supp y) < rat.rk(r) + 1. 
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Proof. By the previous results, y coincides with its first u terms ip(x) = 
Y^=Q c i x ^- Taking a rational 7 < ho and replacing P by 7 P we may 
assume that ho > and that 7 P G T^r* and rat.rk(T*) = rat.rk(r), for 
another finitely generated semigroup T*. 

Define P = P, P i+1 = P;[c^ + Y], T = V* and T i+1 = T { + ^N. 
The coefficients of Pi belong to TZr i - Notice that one has dim<Fj + i> < 
dim <Tj> + 1 and the inequality holds only if Hi G" <rj>. 

For each i, the line L(Pf, Hi) corresponds either to a vertex or to a side 
of jV(Pj). If it corresponds to a side, then there are two different points 
(a, a) and (/?, 6) in C(P«) lying on L(Pj; /Uj). This implies that a, /3 6 Tj and 
Hi = (f3 — a)/(a — b) G <rj>. Hence it is enough to prove that if for an index 
i, [ii corresponds to a vertex of AA(Pj), then for all j > i, [Xj corresponds 
to a side of J\f(Pj). Assume then that Hi corresponds to a vertex v = (a, h) 
lying on L(Pf,Hi)- 

We delay to the following Remark proving that v' = Q(Pi + ±; Hi+i) has 
ordinate equal to one and that ^ , (p i+i;v ')(q ,/i ) 7^ 0, for fi > Hi- 

We can thus write v' = (a' , 1), so that v' is the pivot point of P with 
respect to ij)(x). We know that a' G Tj+i and that, for j > i, 

= 'I'/.,, = V^^te + Coeff(P,;x a '+^ya). 

Since ^ , (p j . t ,/)(g 1 ^) = VE^p.^.^gW) / and Cj / 0, this implies a' + Hj £ r?. 
Hence, fij £ <r«+i> for j > i, which proves that rat.rk({/i; | I > 0}) < 
rat. rk(r i+ i) < rat. rk(T ) + 1. □ 

Remark 4. (Proof that the pivot point of P with respect to ip{x) has ordinate 
1 and that \l/(p .^(g^) 7= 0, for h 7^ Hi)- We write H-, c and P instead of 
Hi, Ci and Pj. Since L(P;h) l~l N(P) = {v}, we have that <E>(p ;/1 )(c) = 

^(p-,v)(Q li )^- Let M(y ,yi) = Ar a y ft + Px a Y ft_1 Yi, g c, be the 

sum of the terms of P corresponding to -u. Then fy(p. v ^(T) = A + BT, 
so that g" = -A/B. Let P = P[c^ + Y] and M =' M[cx^ + Y]. By 
direct computation, M = c h ~ l x a +^ h - 1 ) (AY + BY ± ) + M', where M' has 
only terms of total degree at least two in Yq,Y\, so that the point v' = 
(a + n(h-l),l) is a vertex of M(P) and ^ (py) (T) = c h ' 1 ^ {P . v) (T). Since 

\q\ 7^1, for h'>^ ^ (P y)(^')/0. 



4. g-GEVREY ORDER 

Throughout this section we assume that \q\ > 1. In this case, we prove 
some properties about the growth of the coefficients of a formal power series 
solution of a (/-difference equation. 



Definition 5. A formal power series X^o c * x * ^ s sa ^ ^° ^ e °f Q-Gevrey 
order s > if the series Ylt^o c i M" 5 " x% ^ as a positive radius of conver- 
gence. 
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We will say that a series P = Yla p P<*,p x ° Y p S C[[x, Yq, Y\, . . . , is 
of q-Gevrey order s > if the series 

P*, P \q\-5 s ( a +\ptt 2 x a YP 

(o,p)eNxN n+1 

has a positive radius of convergence at the origin of C n+2 . 

We remark that g-Gervey of order means convergence. This section is 
devoted to proving the following result (the number s(P; y(x)) in the state- 
ment is introduced in Definition [6] and can be computed from the relative 
pivot points of P with respect to y{x)). 

Theorem 4. Let P € C[[x, Yq, Y±, . . . , Y n ]] be a non-zero formal power series 
of q-Gevrey order t > and y(x) £ C[[x]] a solution of P[y] = 0. Then y(x) 
is of q-Gevrey order t + s(P; y) (see the following definition). 

Definition 6. Let Q = (a,h) be the pivot point of P with respect to y(x). 
The number s(P; y) is defined as follows: 

Case h = 1. Let Qj = (aj,hj) be the pivot point of P with respect to 
y(x) relative to the variable Yj (for < j < n). Since Q has ordinate 1, 
Q = Qj f or some j. Let r = max{j | Qj = Q}. There are three cases: 

(RS-R) Ifr = n, then s(P;y(x)) = 0. 

(RS-N) If r < n and hj > 1 for all r < j < n, then s(P;y(x)) can be taken 
as any positive number. 
(IS) If r < n and hj = 1 for some r < j < n, then s(P;y(x)) = 

Case h > 1. By Proposition^ there exist derivatives H = d gyf > w ^ n 
\p\ = h — 1, such that the pivot point of H with respect to y(x) is equal to 1. 
Define s(P;y(x)) as the minimum of all those s(H;y(x))). 

Remark 5. When Q has ordinate h > 1, the number s(P;y(x)) can be 
described directly in terms of the relative pivot points: let Q = (a, h) be 
the (general) pivot point of P with respect to y(x), and Q p (P;y(x)) = 
(a p ,h p ). Let A be the set formed by those 3-tuples (p,i,j) satisfying the 
following properties: \p\ = h, Q p = Q, < i < j < n, and h p > = h, 
where p' = p — &{ + ej. If the set A is empty, we define s(P; y(x)) as any 
positive real number. Otherwise, s(P;y(x)) is the minimum of ^~l a , for 
those (p, i,j) E A. 

Remark 6. Zhang's paper [33J deals with the case in which P is a convergent 
series. The bound given there for the g-Gevrey order of the solution coincides 
with the one described here in cases (RS-R) and (IS), provided h n = 1. In 
the other cases, Zhang proves that some bound exits but without a detailed 
control. In particular, our bound in case (RS-N) is more accurate because 
we prove that the solution is of g-Gevrey order s, for any s > 0. If h n = 1, 
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the bound found in [33J is described with the aid of the Newton-Adams 
Polygon (see [HE]) of the linearized operator along y(x): 

L y = itw [yix)]aJ eC[[x]][(7] - 

3=0 3 

By Proposition [TJ we know that L y is not identically zero if and only if 
the pivot point of P with respect to y(x) has ordinate one. The Newton- 
Adams Polygon N q (L y ) of L y is defined as follows: for each < j < n, 
let lj = ord G NU {oo}. Notice that lj = dj if hj = 1. Then 

Mq(Ly) is the convex hull of the set {(j,lj + r) \ lj ^ oo, r > 0}. It is easy 
to check that s(P;y(x)) is the reciprocal of the minimum of the positives 
slopes of M q {L y ). 

Remark 7. The labels (RS-*) and (IS-*) in Definition [6] correspond to the 
singularity type of the linearized operator L y (regular or irregular). The 
labels (*-R) or (*-N) denote whether the solution y(x) is a regular solution 
of P (i.e. h n = 1) or not. 

4.1. Reduction to solved form. In order to prove Theorem [4], we first 
show (in the paragraphs below) that we may assume that the equation 
P[y] = is in solved form and that the general and all the relative pivot 
points with respect to the variables Yj are reached at step 0. 

Let y(x) = c * x% ^ be a solution of P[y] = 0. We apply the 

process described in the proof of Lemma [7] to P and y(x) in three steps: 

(a) Replace P by some of its derivatives H such that the ordinate of the 
pivot point of H with respect to y(x) is equal to 1 and s(P; y{x)) = 
s(H;y(x)). 

(b) Let be large enough so that all the relative points Qj(H;y(x)), 
for < j < n, have been reached at step N — 1. 

(c) Let 7 = N- 1 and consider P* = W N and y*(x) = J2Zn c * x^ n+1 . 
Then, P*[y] = is in quasi-solved form and P*[y*(x)] = 0. 

If y(x) = YH?=n c i x% i then the relative pivot points of y{x) with re- 
spect to H are the same as the relative pivot points of y(x) with respect 
to H^. Hence, s(H;y(x)) = s(Hj^;y(x)). Finally, by the properties de- 



scribed in subsection 3.1.1 r(Qj(H^;y(x))) = Qj(P*;y*(x)), where r is 
a plane affine map whose restriction to the line of points with ordinate 
one is a translation, so that s(H^;y(x)) = s(P*;y*(x)). This proves that 
s(P;y(x)) = s(P*;y*(x)). Moreover, the general and relative pivot points 
Qj(P*;y*(x)) are reached at step 0. It is straightforward to prove that if 
P is of g-Gevrey order t, then H, Hn and P* are all of g-Gevrey order t. 
Also y*{x) and y{x) have the same q-Gevrey order. This shows that it is 
enough to prove Theorem [4] when the ^-difference equation P[y] = is in 
quasi-solved form and the relative pivot points Qj(P;y(x)) are reached at 
step 0. 
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Finally, assuming that P[y) = is in quasi-solved form, since \q\ > 1, 
the set S + is finite. Let N be an integer greater than the maximum of E + , 
P* = N (P N+1 ) and y*(x) = Yi=N+i d xi ~ N ■ li is clear that s(P;y(x)) = 
s(P*;y*(x)), and also that P* and y*(x) are of the same g-Gevrey order as 
P and y(x) respectively. From this we conclude that we may assume the 
^-difference equation P[y] = is in solved form. 

4.2. Recursive formula for the coefficients. Let y(x) = YaLi c « x% be a 
power series solution of the g-difference equation P[y] = 0, where 

P= E P a>P x a Y''€C[[x,Ya,Y 1 ,...,Y n ]]. 

(a,p)6NxN"+ 1 

Assume that it is in solved form and that the general pivot point Q with 
respect to y(x) and the relative ones Qj = (aj, hj) are all reached at step 0. 
Since the equation is in solved form, Q = (0, 1). Let r be the maximum index 
j i < j < n, such that Qj = Q and let ^(T) be the indicial polynomial of 
P at point Q. From equation (11) one has 

(13) *(g*)ci = -Coeff (Pf, x { Y^), i>l. 

As usual Pi = P[c\x + • • • + Ci-\x l ~ l + Y}. We are interested in computing 
Coeff (P^ x l Y-) in terms of ci, eg, . . . , To this end, we shall consider 
formal series H l in the variables T a p , Cji, x and Yq, Y%, . . . , Y n , where a G N, 
p = (po, . . . , p n ) G N n+1 , < j < n and 1 < I < i — 1, defined as follows 

/i-l \ Pi 

h i = e r «**" n e^ 

(a,p)eN"+ 2 0<i<n \l=l 

For (/3, 7 ) G N x N n+1 , let be the coefficient of a^Y"? in ff\ It is a 
polynomial with coefficients in N and in the variables T a ^ p and Cj i. Denote 
Li = H^q, i.e. the coefficient of x % Y- in IP. A simple computation shows 
that 

L i= E B a,P,d T a,P Y\ IT C U ' 
(a,p,d)ef, 0<j<iJl<Ki-l 

where -B^ d are non-negative integers and the summantion set Ti comprises 

those (a,p,d), such that a G N, p G N n+1 , d = (d i>{ ) G N^ 1 ^" 1 ), for 
< j < n, 1 < I < i — 1, for which the following formulas hold: 

(14) a + E^j,i = *i 

(15) E^ = and so ' E d ^ = i'°i- 

Remark 8. Notice that, substituting in H 1 the variables T QjP for P ajP and 
Cjy for Cjj := q* 1 c\, one obtains Pi. Hence, 

Coeff (P i; ^yQ) =Li(P a>p ,c jtl ). 
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However, in order to have an optimal control on the g-Gevrey growth, we 
need to be more precise and use the position of the relative pivot points of 
P with respect to y(x), and refine the summation set: let T[ be the subset 
of T{ composed by those (a,p,d) satisfying the following properties: 

(16) If j > r, hj > 2, and I > i/2 then d 3 y = 0. 

(17) If j > r, hj = 1, and I > i — aj then dji = 0. 

Let T'l be the complement of J-[ in J-% and let L\ (resp. L") be the sum 
of those terms in Lj corresponding to those (a,p,d) in T[ (resp. in J 7 ','); 
obviously Li = L' i + L'[. 

Lemma 9. The following equality holds: Coeff (P,; x l Y-) = L^P ap ,Cji). 

Proof. Take l with 1 < l < i - 1, and consider P h = P[Y!l=i cj x l + Y] 
(see p])). Let Pi be the series obtained substituting in P\ the expression 



Ya=i x l + Yj for the variable Yj, < j < n. By construction, Pi = Pi, 
hence L { = Coeff (Pi ; x*Y^). 

Write Pi = ^2{ atP )(z Nx ^+^(Plo)a,pX a Y p . Expanding Pi as a series in 
the variables Cjj, lo < I < i — 1, x and Yj, < j < n, let us denote, for 
r < jo < n, by Aj j the sum of terms of Pi in which the variable Cj 0: i 
appears effectively. In order to compute A/ ,z , it is only necessary to take 
into account the terms of P\ in which the variable Yj Q appears effectively, 
that is, only consider the sum over the indices (a,p) G Cj (Pi ). Since we 
are assuming that the pivot point Qj = (aj , hj ) of P with respect to y(x) 
relative to the variable Yj Q is reached at step 0, we may assume that the 
order in x of Aj Q ^ Q is greater than or equal to a JO + hj lo. If jo, lo satisfy the 
premise of either (16) or (17), then aj + hj lo > i and the variable Cj j 
does not appear effectively in the coefficient of x l Y- in Pi . From this one 
infers that L'-(P atP , Cjy) = 0. □ 

From the definition of r, one has \^(T) = Po,e +-Fb,ei T+' ■ -+Po,e r T r , with 
Po,e r 7^ 0. In particular, ^(T) has degree r. Moreover, since the equation 
P[y] = is in solved form, ^(q 1 ) ^ 0, for i > 1. From equation (13) and 
Lemma [9j the following recursive formula holds for all i > 1: 

—1 

( 18 ) a = ^^ L 'i( p ^ c j,i)^ 

where c,y = (f l c\, 1 < I < i — 1 and < j < n. 

4.3. A majorant series. Assume the hypotheses and notations of the pre- 
vious sub-section and that P has g-Gevrey order t > 0. Let s = s(P; y{x)). 
Consider the equation in two variables x and w. 

(19) iu = | g | _a t* |ci|x+ ^ G a , p x a u; |p| , 

(a,p)eC 

where G aiP = |P aiP | | g |-|(«+lpl) 2 +fci(«+lpl)+fc2 ; kl an d ^ 2 are pos itive con- 
stants to be specified later, and C is the set N x N n+1 without the points 



NEWTON-PUISEUX POLYGON OF NON-LINEAR g— DIFFERENCE EQUATIONS 21 

(0,0), (1,0) and (0, ej) for < j < n. It is straightforward to prove that the 
right hand side of (19) is a convergent series and that the equation has a 
unique power series solution w(x) = c 'i x% -> whose coefficients satisfy 
the recursive formulae: 

ci = \q\ r |ci|, c* = Li{G a)P ){c jtl }), i>2, 

where c'- 1 = c[, for 1 < I < i — 1, and < j < n. In particular, c\ > 0, for all 
% > 1, since the coefficients of Lj are non- negative. By Puiseux's theorem, 
the series w{x) is convergent. The following lemma finishes the proof of 
Theorem HI 

Lemma 10. With the above notations, there exist positive constants k% and 
hi such that the coefficients c^ of the solution of equation (19) satisfy 

, s+t ,2 , ,. 

(20) \ci\ < \q\^~ l \c' t \, I > 1. 

Proof. The above inequality holds trivially for 1 = 1. Assume that it holds 
for I = 1, 2, . . . , i — 1. Using equation ( 18 ) and the fact that the coefficients 
of Li are non-negative, one gets 

( 21 ) ^\^y\ E K, P ,,Ga, P , 'ff II (M* M^'Vd ' 



J\d 



■3,1 



= R i( a iPid)B l apA G a , p Y[\ci\ 

(a,p,d)eTi 3,1 

where the indices j and I are < j < n and 1 < I < i — 1, and 

n(a,/9,c0 = / + —l 2 )d u + - (a + \p\f - k x (a + |p|) - k 2 . 
hi 

Claim (proved below): there exist positive constants k\ and k 2 , such that 



(22) Ri(a,p,d) < \q\^\ (a,p,d) € T[ 



Assuming the claim and using equations ( |21[ ) and (22), one gets 

Since the coefficients of Li, the elements G a ,p and are all non- negative 
real numbers, then L'l(G a , P ; {c^}) > 0. Hence, 

^(G Q , p ;{|cl|}) < L<(G a ,„;{c(}) + L?(G a , p ;{c! l }) = Li(G a , p ; {cj}) = K|, 

which proves the Lemma. □ 
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Proof of Claim. Since the degree of \£(T) is r, \q\ > 1 and ^f(q l ) 7^ for 
i > 1, there exists a constant i'fi > 0, such that \q\ tr < K\ | ^(q*) | , for all 
i > 1. Thus, it is enough to prove that there exist k\ > and &2 > such 
that ri(a, p, d) < ^i 2 + ri, for alii > 1 and all (a, p, d) G T[. Grouping the 
terms of and rearranging, we divide the inequality above into two parts 
so that it is enough to prove the existence of positive constants k\ and k 2 , 
such that for all (ct,p,d) £ T[ and i > 1, the following inequalities hold: 



(23) 



(24) 



3,1 3,1 

t , „Q t 



|i 2 + ri + k 2 , 



l 2 d jt i + -(a + |p|) 2 < -* 2 + fci(a + \p\) + fc 2 . 



We first prove the existence for (23) and then for (24). 



Proof of inequality (23). Call r^(a,/?,cQ the left hand side of (23). Let 
T[ = F\ U F2, where iq is the subset formed by those (a,p,d) such that 
I > i/2 implies djj = 0, and F2 is its complement in T[. We shall bound r\ 



in each of F±, F 2 by a polynomial f'(i) = f' 2 i 2 + f[i + f , such that, either 
f 2 < § or f 2 = I and < r. Adjusting k 2 conveniently, one gets (23). 

Let (a,p,d) E F±. This implies that if dji 7^ 0, then I < i/2. As j < n, 
and ^ • ^ / djy < i (which follows form ( |14[ ) ) , we conclude that 



n 1 



r'{%). 



r 'i = 2^ l2dj ' 1 + d M - T ^ + n ^ 1 dj ' 1 ~ ~/ + 

3,1 3,1 3,1 3,1 

If s 7^ 0, then f 2 < s/2. Otherwise, s = 0, and by Definition [6j r = n, hence 
f\ < r. This proves that the polynomial f(i) satisfies our requirements. 
Let (a,p,d) S F 2 . There exists a pair (jo , ^0) such that Iq > i/2 and 



dj 0t i 7^ 1. By inequality (14), this pair is unique and dj j = 1. In this case, 
equation (fl4|) reads as 



(25) a + I + lo = h and in particular / djj, < a, 
j,tylo j,tyio 

where a = i — lo < i/2. This implies also that for I 7^ Iq and dji 7^ one has 
I < a. Therefore, 



Y l2d 3\i + l o \ + Y jl d hi + h k 



d,¥lo 



j,¥io 



< - j a Y 1 d j,i + (i ~ a ) 2 + n Y 1 d i> 1 + ~ a ) 
h¥lo J 3,Wa 



< -^{2a 2 — 2ai + i 2 ) +na + jo(i — a) 

= (s/2)i 2 + (jo - sa)i + (sa 2 + na - a j ) := /»(»). 
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For a fixed i, the graph of fi(a) is an upwards parabola, so its maximum in 
an interval is reached at its endpoints. The available range for a depends 
on jo. If jo < r, then a £ [l,i/2[, and we take f'(i) = s/2i 2 + ri + fg. We 
can chose f' in such a way that max{/j(l), /j(i/2)} < f'(i), for all i > 1, 
because the coefficient of in fi(i/2) is s/4 < s/2 and the coefficient of 



2 in fi(l) is jo — s < r. If, on the other hand, jo > r, case (16) does 
not happen because lo > i/2 so that hj = 1 and lo < i — aj , and the 
range for a is [aj ,i/2[. By definition of s, one has jo — so JO < r and 
s > 0. Consider r'{i) = s/2i 2 + ri + f , where f is chosen in such a 
way that max{/j(?'/2), fi(dj); j > r,hj = 1} < f(i), for all ? > 1. Such 
an f' exists because the coefficient in fi(i/2) of i 2 is s/4 < s/2, fi(aj) < 
(s/2)i 2 + (j — s a,j)i + s a 2 + na^ and j — so, < r for those j such that j > r 
and hj = 1. 



Proof of inequality (24). For t = 0, the inequality holds trivially, so we 
may assume that t > 0. Let (a,p,d) G J 7 ^. Denote = ^2™ =Q dj : i, for 
1 < Z < i — 1 and let lo be the maximum of the indices / such that di ^ 0. 
From equations ( [T4| ) and ( [To] ), the fact that Iq > 1 and di > 1, one gets: 

i — \p\ = a + / (i/ — = a + y~](/ — l)di + (lo — 1 ) <i; > a + lo — 1. 



From which /o > a+|/?| — 1. Taking into account that a > 0, equation (14), 
and the fact that lo > I for any I with di ^ 0, we conclude that 

i 2 = (i-lo + lo) 2 = {i-lo? + ll + 2lo (i-l ) 
> (a + \p\ -l) 2 + l 2 o + 2l lj2 ld i + l o{d lo -l) 

>(a + \p\ -l) 2 + l 2 + Y^ l 2 d t + l 2 (d lo - 1) 

Wo 

>{a + \p\) 2 -2{a + \p\) + Y,l 2 d l . 

i 

This gives inequality ( 24 ) for k\ > 2/t and finishes the proof of Theorem |4j 

□ 



5. Working example 

Let us consider the g-difference equation P[y] =0 of order 5 and degree 
6, where 

P = 4Y l *-9Yo 2 Y l Y 2 + 2Yo*Y 2 -x*Yo 4 Y^ + 3 ^-^-x*Yo + x 5 , 

and q = 4. Its Newton Polygon is Af(P) in Figure [TJ It has four vertices 
vo = (3,6),t>i = (0,4),f2 = (1)2), «3 = (5,0) and three sides L\,L 2 and L3 
with respective co-slopes 71 = —3/2, 72 = 1/2 and 73 = 2. We apply some 
steps of Procedure [l] to P. As P is a polynomial, = —00. 
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Figure 1. Newton Polygons from the working example. 



In order to find all the possible starting terms cq of a solution, we 
need to consider all the vertices and sides of Af(P) according as formulas ^ 
and (g). For the vertices, we get: * (P . Vo) (T) = -3T 10 , ^ (P . Vl) (T) = T 2 (T- 
2)(4T- 1), ^(p ;V2) (T) = T 2 /? 4 , * ( p^ 3 )(T) = 1. Hence, for j = 0,1,2,3 the 
only satisfiable formula in ^ is the one corresponding to vertex v\, that is 
^ (P;vi)(Q ,i ) = an d —3/2 < fx < 1/2. This gives /i = — 1 for any nonzero 
c. For the sides, we get: $(p ;71 )(c) = q~ lb c 4 {2q 12 - 9q 21/2 + 4q 9 - c 2 ), 
$(P ;72 )(c) = c 2 /64, and <&(p ;73 )(c) = (c — l) 2 . According as ([5]), the only 
possible starting terms related to the sides are ±1024\/l5 x~^l 2 and x 2 . 
Notice that L2 gives rise to no starting term. 

Following Procedure [l] we choose x 2 , that is Co = 1 and (jlq = 2. The 
polynomial Pi = P[x 2 + Y] has 33 terms that we do not exhibit; its Newton 
Polygon is J\f(Px) in Figure [Tj Since y = is not a solution of -Pi [2/] = 
because C(Pi) has points on the OX-axis, we need to perform step (a. 2) of 
Procedure [TJ that is finding /j, > /xo = 2 and c / so that ^(p i;M )(c) = 0. 
Thus, we can only use the vertices V2 and v' 3 and side L' 3 . 

For formula ^ we get ^(p i; „ 2 )(T) = ^^ P . V2 ^ and that ^^.^(T) is a 
constant, hence those vertices do not give rise to subsequent terms. For side 
L' 3 , we get [i\ = 7/2 and ^/(p i;/J1 )(c) = 64c 2 + 225792, so that there are two 
possibilities for c\. We choose c\ = 21\/8\/— T and go on with Procedure hi 
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Let us consider P2 = Pifcix^ 1 + Y] whose Newton Polygon is N(P2) 
having a side L 3 ' of the same co-slope as L' 3 and another L4 of co-slope 5. 
As vertex gives ^ r (p 2;t ,//)(g /i ) = q 2 ^ + 16384 which has no real solutions, it 
is useless to find fi2- Hence we must use L'l which gives \xi = 5 and (after a 
trivial computation) C2 = —88984/65. 

Notice that, after performing the first two steps detailed above and getting 
x 2 + 21 v / 8\/ r T2; 7//2 , the fact that Vg gives rise to a formula which no fi > 7/2 
can satisfy and that it has ordinate one implies that, taking P* = ^^2, the 
equation P* [y] = is solved form. Therefore, by Lemma [8] there exists a 
unique solution of P[y] = of the form: 

y(x) = x 2 + 21v / 8v /Z T x 7/2 + o(x 7/2 ). 

Notice also that as P* G Cffo; 1 / 2 ]] [Y], Lemma [8] guarantees as well that 
y[x) G C[[x 1 / 2 ]}. 

The pivot point of P with respect to y(x) is Q(y(x); P) = v% = (4.5, 1) and 
Y2 is the highest order appearing effectively in it, hence r = 2 in Definition [5] 
There being no monomials with 13 or 14 in P we only need consider the 
pivot point relative to I5 which is the point Q e5 (y(x); P) = (13, 1) (notice 
that C es (P2) is in the region above and to the right of the dashed line). 
Applying Definition 5] formally we would get s(y(x); P) = ^-45 = 6/17. 

Concerning the growth of the coefficients of y(x), we transform it into 
a formal power series in order to apply Theorem |4j We do this by means 
of the ramification x = t 2 . The series y(t) is a solution of a ^-difference 
equation P[y] = derived from P with q = q 1 ! 2 . The ramification induces 
a horizontal homothecy of ratio 2 on the cloud of points of P, P\ and P2- 
Hence s(y(t); P) = 2(13-4 5) = 3/17 is a bound for the g-Gevrey order of y(t). 
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